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In a two-dimensional two-component plasma, the second moment of the number 
density correlation function has the simple value {127r[l — where T is 

the dimensionless coupling constant. This result is derived directly by using dia- 
grammatic methods. 



Dedicated to Joel Lebowitz on the occasion of his 70th birthday 



1 Introduction 

The system under consideration is the two-dimensional {2d) two-component 
plasma (TCP) (Coulomb gas), i.e. a system of positive and negative point- 
particles of charge ±g, in a plane. Two particles at a distance r from each 
other interact through the 2d Coulomb interaction =i=g^ ln(r/L), where L is 
some irrelevant length (set to unity, for simplicity). Classical equilibrium sta- 
tistical mechanics is used. The dimensionless coupling constant is F = jSq"^, 
where /5 = 1/ kT is the inverse temperature. The system is known to be stable 
against collapse when F < 2. Let {n{0)n{r))T be the truncated density-density 
correlation function, where n{r) is the total microscopic number density at r. 
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In a recent paper [1], an exact expression for the second moment of that 
density-density correlation function has been derived: 



However, the two ahernative derivations which were given were rather indi- 
rect ones, and the Conclusion of the paper had the sentence "A more direct 
derivation is still wanted" . At about the same time, another paper [2] used a 
diagrammatic method for deriving an exact expression for the sixth moment 
of the pair correlation function of the 2d one-component plasma (OCP). The 
authors of both papers soon realized that the diagrammatic method of Ref. 2 
could be easily adapted to the two-component case and used for providing a 
more direct derivation of (1), which is the main content of the present paper. 

The paper is organized as follows. In Section 2, we generalize a renormalized 
Mayer expansion [2], [3], to the case of a many-component fluid. The appli- 
cation of the formalism to the 2d TCP in Section 3 leads, with the aid of the 
results obtained in Ref. [2], to the sum rule (1) for the density correlation func- 
tion. The charge correlation function is discussed in Section 4. Section 5 deals 
with consequences on finite-size effects for the sphere geometry, confirming 
predictions from conformal invariance [4], [5]. 



2 Renormalized Mayer expansion for many-component fluids 

We first consider a general classical multi-component fluid in thermody- 
namic equilibrium. The last system is composed of distinct species of par- 
ticles {a} with the corresponding densities {n{r,a)} at vector position r 
of a d-dimensional space. The particles interact through the pair potential 
v{ri,ai\rj,aj) which depends on the mutual distance of particles i,j as well 
as on their types cxj, aj. Vector rj will be represented simply by i. 

According to the Mayer diagrammatic formalism in density [6], the gener- 
alization from one-component to more-than-one component fluid consists in 
considering an additional particle-type label a for each, fleld (black) or root 
(white), vertex. Simultaneously, besides the integration over spatial coordi- 
nate of a fleld circle, the summation over all cr-states at this vertex is assumed 
as well. The renormalized Mayer representation of the excess Helmholtz free 
energy [2], resulting from the expansion of each Mayer function in the in- 
verse temperature and the consequent series elimination of two-coordinated 
fleld circles between every couple of three- or more-coordinated fleld circles 
(by coordination of a circle we mean its bond-coordination, i.e. the number of 
bonds meeting at this circle), can therefore be straightforwardly extended to 
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the many-component case. The renormahzed i^-bonds are now given by 
K -(3v 

0^y^,N^^O = O- -O + O • -O + O- • • -O + ... 

l,ai 2,(72 2,(72 1,0-1 2,(72 1,0-1 2,(72 

or, algebraically, 

X(l,ai|2, (72) = [-/5t;(l,ai 12,^2)] 

+ E /Mt^(l,ai|3,a3)] n(3,a3) K(3, a3|2, ^2) d3. (2) 

The procedure of the bond-renormalization transforms the ordinary graphical 
Mayer representation of the free energy into 

00 

(3F'^[n] = H, + ^ D^'\n\, (3a) 

3 = 1 

where is minus the excess free energy, 

A f- - - 

D(o)= / \ + ; ; + ... (3b) 

*^ — ^ 4----* 

is the sum of all unrenormalized ring diagrams (which cannot undertake the 
renormalization procedure because of the absence of three- or more-coordinated 
field points) and 

00 

^ = |all connected diagrams which consist of iV > 2 field 

s=l 

n(i, (7j)— circles of (bond) coordination > 3 and multiple 

K {i,ai\j, a j)— bonds, and are free of connecting circles | (3c) 

represents the set of all remaining completely renormahzed graphs. By multi- 
ple i^-bonds we mean the possibility of an arbitrary number of i^-bonds be- 
tween a couple of field circles, with the topological factor 1/ (number of bonds)!. 
The order of s-enumeration is irrelevant, let us say 




JJ(^) ^(4) ^(5) 



etc. 

The truncated pair correlation 



l^f, 10 N ^2(1,^1 2,(72) T 

h{l, cxi\2, (X2) = — -— r - 1 (5) 

n{l,ai)n{2,a2) 



is related to the pair direct correlation c via the Ornstein-Zernike (OZ) equa- 
tion 



/i(l,ai|2,a2) = c(l,o-i|2,cr2) 

+ H / c(l, o-i|3, 0-3)72(3, o-3)/i(3,o-3|2,a2)d3. (6) 

The free energy is the generating functional for c in the sense that 

c(l,cri|2,a2) = — (7) 
dn{l, ai)dn{2, a2) 

According to (3), the direct correlation is thus written as 

00 

c(l,ai|2,a2) = o - - o +c(°)(l,ai|2,a2) + $]cW(l,^i|2,f^2), (8a) 

2,(72 s=l 



where c^^\l, ai|2, (T2) = S'^D^'^^ Sn{l, ai)Sn(2, (T2) can be readily shown to cor- 
respond to the renormalized Meeron graph 

cW(l,cri|2,a2)=^<3:j)^ = ^K\l,a,\2,a2) (8b) 



and 



52 D(^) 

c(^)(l,ai|2,a2) = — (8c) 

dn{l, ai)dn{2, a2) 



(s = 1, 2, . . .) denotes the whole family of (1, <Ji){2, cr2)-rooted diagrams gen- 
erated from D^'^\ To get explicitly a direct-correlation family, one has to take 
into account the functional dependence of the dressed i^-bonds (2) on the 
density, too. It is straightforward to verify by variation of (2) that it holds 

6K{l,ai\2,a2) , 

— ^— r — = K{1, ai 3, a3)K{3, ^3 2, a2). (9) 
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As a consequence, the functional derivative of D^'^^ with respect to the den- 
sity generates root-circles, besides the field-circle positions, also on i^-bonds, 
causing their "correct" K-K division. For example, in the case of the generator 
D^^^ drawn in (4), we obtain 



2,(72 1,0-1 




3 Density correlations in the 2d TCP 

Let us now concentrate on the neutral 2d TCP of positive {a = +) and 
negative (a = — ) charges, with the Coulomb interaction energy given by 

- Pv{i,ai\j,aj) = criajT\n\i - j\. (11) 



We consider the infinite- volume limit, characterized by homogeneous densi- 
ties n(l, 0") = Her and both isotropic and translationally invariant two-body 
correlations, c(l, cri|2, (T2) = c^^aadl - 2|), /i(l, cri|2, (T2) = /if^i^^jdl - 2|). The 
requirement of the charge neutrality, n+ = n_ = n/2 (n is the total number 
density of particles), automatically implies the symmetry with respect to the 
-|- ^ — interchange of the charge, 

= (= h+), h+_ = h_+ (= /i_); (12a) 

C++ = c__ (= c+), c+_ = c_+ (= c_). (12b) 



Thus, in Fourier space 

fir) = ^ I expilk ■ r)f{k)d'k, (13a) 



fik) 



27r 



exp(— zfc- r)f(r)d^r 



1 

2^ 



/(r)r'M>, 



(13b) 
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the OZ equation (6) implies 



h_{k) = c_{k) 1 + imh^{k) + nnc^{k)h_{k), 
where k — \k\. Introducing 



(14a) 
(14b) 



a,a'=± 

2^ C^a' = -(C+ + C_), 
a,a'=± 

the summation of (14a) and (14b) yields 

h{k) = c{k) + 27mc{k)h{k). 



(15a) 

(15b) 

(16) 



In the renormalized format, due to the translational invariance of renor- 
malized bonds K{l,ai\2,a2) = K^^^^^il^ ~ 2|), the Fourier transformation of 
relation (2) results in 

Kaa'ik) = [-Pvaa'ik)] + 27r ^ [-/9w'(A:)] n^''K^''a'{k). (17) 

o-"=± 



With regard to 



-(3vaa'(k)^-aa'r/k\ 



(18) 



equation (17) yields 



Kaa'{k) = -aa'- 



r 



k^ + 2nr{n+ + n_) ' 



r J27rr(n+ + n_) 



(19a) 
(19b) 



with Kq being the modified Bessel function of second kind. The pair direct 
correlation (8) is expressible as 



-'ITl(T2 



(|1 - 2|) = (7i(72ln|l - 2| + ^ [ri^o(|l - 2|V27rrn)'' 



+ E 



.is) 



c^7..(|l-2|)- 



(20) 



s=l 
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c, equation (15b), is then given by 



1 n 2 °° 

c(|l-2|) = - ri^o(|l-2|V'2^) +Y.c^'\\l-2\), (21a) 



where 



c(^)(|l-2| 



E 4Kli-2|) 



0-l,(T2 = ± 



_ 1 



n{r,a-)=n/2 



(21b) 



As concerns the evaluation of because of the symmetry Xo-^g. = —K^r-a, 
it is rather obvious that for n+ — n_ the summation on signs gives a non- 
vanishing c^^^ only if all the vertices of the generator D^^^ have an even bond- 
coordination [note that after the functional derivation (21b), with regard to 
(9), all the white circles will also have an even bond-coordination]. Further- 
more the resulting c^^^ is the same as in the case of the OCP with particle 
density + n_ = n. It was proven in Ref. [2] that the zeroth- and second- 
moments of such a c^^^ vanish independently of the topology of D^^\ so that 



(22) 



To summarize, using formula [7] 

oo 



(2s + l)! 



(s>0) 



(23) 



in the Fourier representation (13b) of and regarding equation (22), the 
small-/c expansion of the Fourier component of c (21) reads 



c{k) 



fc2 



87m 96(7rn)^ 



(24) 



Inserting (24) into (16), we find 
h{k) ^ 



k^ 



Sttu (1 - 5) 96(7rn)2 (l - f )' 



(25) 



The truncated density-density correlation 
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is related to h [see definitions (5), (12a) and (15a)] via 

{h{0)h{r))T = n^h{r) + n5{r). 



(26) 



With respect to (13b), we finally obtain 

r 



n[h{r)d'r^ (27a) 
^ 4(1-0 

/ {niO)n{r))T dV= ^. (27b) 

127r (l - 



The first formula (27a) is identical to the compressibility sum rule 

dn 

gjp) 



with the use of the exact equation of state (3p = n(l — r/4), where p is 
the pressure. The second formula (27b) coincides with the new sum rule (1) 
suggested in Ref. [1]. 



4 Charge correlations in the 2d TCP 



One might try to use the present formalism for studying the charge corre- 
lations. Unfortunately, one can only reproduce known results. Indeed, intro- 
ducing the charge correlation and direct correlation functions 

h' = ^{h+-h.), c' = ^(c+-c_), (29) 



one now finds that the non-vanishing graphs contributing to c' arc the ones 
which have their two root-vertices with an odd coordination and their field 
vertices still with an even coordination. For instance, only the first diagram in 
the c^^) family [see formula (10)] survives. The argument of Ref. [2] no longer 
applies, and the corresponding contribution to c'(A;) is 0{1) rather than 0{k^). 
One finds 

c'{k) = -T/e + 0(1) (30) 



and this leads only to the known Stillinger-Lovett sum rules for the zeroth 
and second moments of the charge correlation function h'. 
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5 Consequences 



Cranking backwards the arguments in Ref. [1], it is now possible to derive 
conformal-invariance properties from the 2nd moment (27b). Indeed, within 
the scheme estabhshcd in Ref. [1], the logarithm of the grand-canonical parti- 
tion function of the 2d TCP, formulated on the surface of a 3d sphere of radius 
R, can be written as 

InS = 4rf/3p + /(i?), (31a) 



where, in the limit — > oo, 

^ = - 9'^^ / (^iO)n{r))T dV. (31b) 

Assuming the sum rule (27b), one finds the behaviour, as i? — > oo, 

f{R) ~ In i? + const, (32) 



i.e., for an arbitrary coupling F, the finite-size correction to InS exhibits the 
universal form suggested by using predictions of the conformal-invariance the- 
ory [4] , [5] , with the correct prefactor for the sphere geometry. 

The mapping of the 2d TCP onto the massive Thirring model confirms the 
validity of a sum rule for the latter model (for details, see Section 4 of Ref. 

[I])- 



6 Conclusion 



Unlike previously known moments in Coulomb systems, the new 6th mo- 
ment of the pair correlation function for the one-component plasma, and the 
new 2nd moment of the density (not charge) correlation function for the two- 
component plasma have been exactly obtained only for two- dimensions and 
point-charges (no short-range interactions). A simple basic ingredient of the 
derivations is that these systems have only one length scale provided 
by the density; the quantities of interest depend on the density only trivially 
through that length scale. Nevertheless, technically, showing that higher order 
graphs give to c contributions 0{k'^) (Ref. [2]) is somewhat lengthy. Whether 
a simpler derivation is possible is an open problem. 
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